Effect of Hadronic Rescattering on Elliptic Flow Following Hydrodynamics Model 
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The effect of hadronic rescattering on the elliptic flow has been investigated by the Cooper-Frye 
hadronization model from hydrodynamic evolution following by the afterburner hadronic rescatter- 
ing model for 200 GeV/c Au + Au at 20-40% centrality. It is found that the hadronic rescattering 
can suppress elliptic flow V2 in momentum space especially in lower transverse momentum region. 
In addition, the hadronic rescattering effects on transverse momentum spectra and anisotropy co- 
ordinate space of hadrons are studied. 

PACS numbers: 25.75.-q, 24.10.Nz, 24.10.Pa, 25.75.Ld 



I. INTRODUCTION 



Lattice QCD calculations have predicted a transition 
from quark-gluon matter (QGP) to hadronic matter at 
high temperature or density Q, Q . This transition was 
believed to occur around ten microseconds after Big 
Bang. Practically, one may obtain this piece of infor- 
mation by studying the relativistic heavy ion collisions 
at Brookhaven National Laboratory now. Somepossible 
probes have been proposed for searching for it Q • El- 
liptic flow V2 is one of these probes, from which one may 
understand some early information in relativistic heavy 
ion collision through measuring anisotropy of momen- 
tum space of particles in final state. Some experimen- 
tal results from BNL-RHIC have been published recently 

iiaiiimiia. 

As we know, hydrodynamical model, which is one 
method to study relativistic heavy ion collisions, can re- 
produce i>2 data very well below ~ 2 GeV/c of trans- 
verse momentum (p?) 0> El Ejl [IB EH However, the 
hadronical rescattering effect on elliptic flow after the 
hadronization following by the hydrodynamics, at least 
in our knowledge, has not been investigated in details so 
far. Based on this motivation, we shall focus the influ- 
ence of rescattering in hadronic state on elliptic flow V2 ■ 
It should be pointed that hadronic phase space before 
rescattering is produced directly by an anisotropic mo- 
mentum space of produced hadrons in order to produce 
v 2 in our calculation, which will be described in detail in 
following. As a result, we find that v-i of final particles 
is suppressed due to the rescattering among hadrons to 
a considerable extent. It may indicate that the effect of 
rescattering among hadrons on vi is not ignored in the 
research on 
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II. HYDRODYNAMICS MODEL AND INITIAL 
PHASE SPACE OF HADRONS 

Hydrodynamics of relativistic heavy ion collision [TtI 
E EH is defined by (local) energy-momentum and net 
charge conservation, 

= 0, 9 M Yf = 0, (1) 

where denotes the energy-momentum tensor and Nf 1 
the net four-current of the ith conserved charge. Here 
will only take the net baryon number. We implicitly 
assume that all other charges which are conserved on 
strong-interaction time scale, e.g., strangeness, charm, 
and electric charge, vanish locally. The corresponding 
four-currents are therefore identically zero, and the con- 
servation equations are trivial. 

For ideal fluids, the energy-momentum tensor and net 
baryon current assume the simple form: 

r = ( e+P K/- w ^, jv£ = pbu", (2) 

where e, p, and ps are energy density, pressure, and net 
baryon density in the local rest frame of the fluid, which 
is defined by Ng — (pb,0), respectively. Here g^ v — 
diag(+, -,-,-) is the metric tensor, and — 7(1, the 
four-velocity of the fluid ( v is the three-velocity and 7 
= (1 — V 2 ) -1 / 2 , the Lorentz factor). 

For simplicity, the hydrodynamics that we are using as- 
sumes cylindrically symmetric transverse expansion with 
a longitudinal scaling flow profile, v z = z/t though we are 
investigating the non-central collision. We recognize that 
our initial phase space of hadrons in this way is not per- 
fect since we put an initial anisotropy of phase space by 
hand just after Cooper-Frye hadronization as illustrated 
in the following. Nevertheless, we think that it is still 
useful and interesting to investigate the hadronic rescat- 
tering effect on spectra and flow based on this treatment. 
In this point, our simulation is somehow like that we start 
from an assumed initial phase space of hadrons and then 
study the hadronic rescattering effect. This phase space 
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is similar to a sudden hadronization phase space which 
the partonic anisotropy turns to hadron anisotropy. For 
hydrodynamics, at z = 0, Eqs (1) reduces to 

d t E + d T [{E + p)v T ] =-(— + b(E + p), (3) 

Tt t 



Vt i 

d t M + dr[Mv T + p] = -(— + -)M, (4) 

rj 1 t 





d t R + d T [Rv T ] = -( 



)R, 



(5) 



where we defined E = T 00 , M = T 0T , and R = N%. In 
the above expressions, the index T refers to the transverse 
component of the corresponding quantity. The equations 
(3)-(5) describe the evolution in the z = plane. Because 
of the assumption of longitudinal scaling, the solution at 
any other z / can be simply obtained by a Lorentz 
boost. 

Initial conditions for Au + Au at 200 GeV/c in 20- 
40% centrality is supposed that initial energy density = 
10.3 GeV/ fm? and initial time = 0.6 fm/c, and a hy- 
drodynamical evolution of cylindrically symmetric trans- 
verse expansion with a longitudinal scaling flow profile is 
adopted. When temperature of liquid reaches to a criti- 
cal energy density ( = 0.42 GeV/ fm 3 ), a hadronization 
of QGP which is described by employing the MIT bag 
model EOS [2lJ will take place. And an ideal hadron 
gas is established and the corresponding EOS is ap- 
plied. Hadronization in hydrodynamics usually takes the 
method of Cooper and Frye 22] . In this work, however, 
we make an assumption that the coordinate and momen- 
tum spaces of hadronized particles are not isotropic but 
have shapes as plotted in Figure 1 in order to simulate 
Au + Au at 200 GeV/c in 20-40% centrality and get v 2 
in our calculation. Here the mean anisotropy in momen- 
tum space, F =<Px/Py>, takes 1.18 in our simulation 
for all produced hadrons. 

For the details of the simulation to result in the above 
anisotropic distributions, we wrote 

R x = Rt x cos9, P x = Ptx cos?/> 



R y = RTySinO, P y = Pr y simp 



R z = t * tanhy, P z = mx sinhy 



FIG. 1: Contour plots of shapes of the phase space where hadrons 
are produced when two same spherical overlap in half-central col- 
lision. Left: transverse plane at z = in coordinate space; Right: 
transverse plane at Pz = in momentum space. 



< 0/v >= 1 



(8) 



where tut = \/Vt + TO 0' (Rx,Ry,Rz,t) and 
(P x , P y , P Z ,E) are the 4- vectors of produced hadrons 
in coordinate and momentum spaces respectively. We 
assume that hadrons are randomly distributed in the 
overlap coordinate space of two spherical nucleus in 
terms of the left column of equations (7) by hadroniza- 
tion. £ and y are random numbers which are respectively 
randomly distributed in [0,14] fm/c and [-5,5]. As for 
momentum and position of hadrons, Pt x and Rt v denote 
a random variable distributed in [0,4] GeV and [0,7] fm, 
then Pt v and Rtx can be obtained by the anisotropic 
factor in momentum and coordinate space, namely 
F, which makes our momentum and coordinate space 
distribution in transverse plane elliptic and anti-elliptic. 
If we make F equal to 1, our hadronization method will 
go back to the method of Cooper and Frye. 9 and ip 
are the azimuths distributed in [0,27r] in coordinate and 
momentum spaces. But here we assume these produced 
hadrons are from a radial source in which and "0 have 
a correlation like the equation (8), which means that 
once a tj) has been decided, the corresponding 9 will 
be selected in a Gaussian distribution whose mean is 
ip and whose width is always assumed to be 0.27T in 
our calculation. The set of equations results in the 
coordinate and momentum space as figure 1 shows. 
Here we choose the anisotropic factor F in momentum 
space, equals to 1.18 for 200 GeV/c Au + Au in 20-40% 
centrality. 



t = dtanhy, E = yj P* + P* + P* + m* (6) 
< Ptx/Ptv >=< Rtv/Rtx >= F (7) 



III. HADRONIC RESCATTERING MODEL 

After the hadrons are produced, they enter next rescat- 
tcring process among them. We here use the rescattering 
model from LUCIAE model HHH]. Two particles will 
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collide if their minimum distance d mm satisfies 



&total 



TV 



(9) 



where atotai is the total cross section in fm 2 and the min- 
imum distance is calculated in the C.M.S. frame of the 
two colliding particles. If these two particles are moving 
towards each other at the time when both of them are 
formed, the minimum distance is defined as the distance 
perpendicular to the momentum of both particles. If the 
two particles are moving back-to-back, the minimum dis- 
tance is defined as the distance at the moment when both 
of them are formed. Assuming that the hadrons move 
along straight-line classical trajectories between two con- 
secutive collisions it is possible to calculate the collision 
time when two hadrons reach their minimum distance 
and order all the possible collision pairs according to the 
collision time sequence. If the total and the elastic cross 
section satisfies 



^elastic 
&total 



> 



(10) 
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FIG. 2: The transverse momentum distributions of n~ , k~ , P in 
Au + Au at 200GeV/c in 20-40% centrality in our hydrodynamical 
+ rescattcring model. Full lines and dash lines represent pt spectra 
after and before hadronic rescattering, r esp ectively; solid symbols 
are experimental data which come from l2, r ^ . 



where r\ is a random number then the particles will be 
elastically scattered or else the collision will be consid- 
ered as an inelastic collision. The distribution of the 
momentum transfer, t, is taken as 



IV. 



RESULTS AND DISCUSSIONS 



da 

- ~ exp{Bt) 



(11) 



where B, for an elastic scattering, depends on the masses 
of two scattering particles. The azimuthal angle will be 
isotropically distributed. 

During the rescattering process, the following inelastic 
reactions are included in our calculation: 



7rN ^ A 7T 


7rN ^pN 


NN^AN 


7T 7T ==^k k 


7T N ^ k Y 


nN^kY 


7T Y ^ k S 


nY^kZ 


k N ^ 7T Y 


kN ^ ttF 


k Y 7T S 


kF^TT H 


k N v± k 3 


kN ^kE 


7T E ^ k fi- 


7T S k n~ 


kS ^ 7T fF 


k S ^ 7r Q~ 


N N annihilation 


Y N annihilation 



where the hyperons Y = A or E. The relative probabili- 
ties for the different channels, e.g. in (n N)-scattering, is 
used to determine the outcome of the inelastic encounter. 
As the reactions introduced above do not make up the 
full inelastic cross section, the remainder is again treated 
as elastic encounters. 



Figure 2 shows the transverse momentum distribution 
of 7r~, fc~, P before and after rescattering in Au + Au 
at 200GeV/c in 20-40% centrality in our hydrodynamical 
calculation followed by the rescattering model. 

At the same time, we obtain the v 2 of hadrons before 
and after rescattering in terms of formula (12), which 
indicates the anisotropy of hadrons' momentum space. 



1>2 =< 



Px-Pv 



P 2 y 



> 



(12) 



Figure 3 shows dependence of vi of 7r~+fc~, h~+h + , P, 
A+A on pt before and after rescattering in hydrodynam- 
ics and experimental data in Au+Au at 200GeV/c in 20- 
40% centrality p5j . 

In order to study anisotropy of hadrons' coordinate 
space, the parameter named as £2 is defined by formula 
(13). 



£2 =< 



> 



(13) 



The plots (a) (c) and (e) in figure 4 give px spectra, 
V2 and 62 of four types of particles (( ir + +ir~), (fc + +fc _ ), 
(P+P) and (A+A) ) vspr before and after rescattering in 
hydrodynamics+rescattering model. The open and filled 
points, respectively, stand for the status before and after 
hadrons' rescattering. We can see that pr yield and £2 
of hadrons is enhanced by rescattering, and also note £2 
changes its sign through rescattering. On the other hand, 
elliptic flow V2 of hadrons is suppressed by the rescatter- 
ing. In order to quantify the extent of enhancement (or 
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FIG. 3: The px dependences of V2 of ir~+k~, h~ , P and A 
+ A before and after rescattering in hydrodynamical+rescattering 
model and experimental data which come from |3 and ll ll in 
Au+Au at 200 GeV/c in 20-40% centrality. 
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FIG. 4: (a),(c) and (e): The pt dependences of spectra, V2 and £2 
ofn~*~+n~ (square), k + +k~ (circle), P+P (up-triangle) and A+A 
(down-triangle) before (open symbols) and after (full symbols) 
rescattering in hydrodynamical + rescattering model in Au+Au 
at 200 GeV/c in 20-40% centrality; (b),(d) and(f): The p T depen- 
dences of sp T , s V2 and D e2 . See more details in text. 



suppression ) at different px bins, we define three factors 
named as sp T1 s V2 and D e2 by formula (14), (15) and (16). 



dydP T a t ter 

s pt(pt) = str- 

dydP T before 



s V2 (pt) 



V2(PT)afte 
V2(pT)befo'< 



D e 2 (PT) = e(pr) 



after 



T /before 



(14) 



(15) 



(16) 



The plots (b) (d) and (f) in figure 4 show respective 
dependence of sp T s V2 and D e2 onpT- In (b), we find that 
the yield of hadrons is increased by hadrons' rescattering, 
especially for higher pt yield. It indicates the rescatter- 
ing produced many secondary hadrons, whose px seems 
harder than those hadrons in last generation. In (d), s V2 
increases with px, which indicates that the suppression 
effect on vi from hadrons' rescattering becomes weaker 
and weaker with the increasing of px- In (e) and (f), we 
see the rescattering turns the shape of coordinate space, 
and the effect becomes stronger with the increasing of px ■ 
It indicates the rescattering makes hadrons' momentum 
space less anisotropy and more heavier particles, more 
suppression the vi . While the changes of the tropism of 
coordinate space occurs in the same time. In our calcu- 
lation, at least 20 ~ 40% V2 will disappear through final 
hadrons' rescattering. So the effect from final hadrons' 
rescattering on V2 maybe should be considered in order 
to research early state information in relativistic heavy 
ion collisions. 

In order to investigate the initial geometrical depen- 
dence of our results farther, we next attempt to take 
Cooper-Frye hadronization method but with three dif- 
ferent anisotropic conditions to see the change of effect 
from final hadrons' rescattering. Since the status of mo- 
mentum and coordinate space before the hadrons' rescat- 
tering can not be measured by us, we assume anther three 
conditions here. 

(I) We decrease the correlation between momentum 
space and coordinate space, i.e. increase <j= 0.2ir to a 
— 0.5-7T. As we see the circles and squares in figure 5, 
the stronger correlation seems to produce more hadrons 
and suppress vi more strongly, but its change for the 
anisotropy of coordinate space €2 is slight. 

(II) We only remain anisotropy of coordinate space 
when hadronization takes place. What will happen af- 
ter rescattering among hadrons? The diamons in fig- 
ure 5 show our results. Obviously, only hadronic rescat- 
tcrring can not induce V2 if the interaction system has 
not anisotropy of momentum space before the hadronic 
rescattering, though its high sp T indicates that more sec- 
ondary hadrons are produced in rescattering. On the 
other hand, the anisotropy of coordinate space also dis- 
appear after hadrons' rescattering. It is also reasonable 
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FIG. 5: (a), (c) and (e): The pt dependences of spectra, V2 and €2 
of charged particle before (open symbols) and after (full symbols) 
rcscattcring in hydrodynamical + rescattering model in Au + Au 
at 200 GeV/c in 20-40% centrality. The circles and squares are 
for the hadronization with both initial 1)2 and 62, and the circles 
correspond to 0.27T of the width of i/>-<f> correlation and squares 
correspond to 0.57T. The diamonds are from the hadronization 
whose particles are only with initial t2and the triangles are from 
the hadronization whose particles are only with initial V2 (see the 
insert of (b)); (b),(d) and (f): The px dependences of sp T , s V2 and 
D 62 in these four hadronaztion situations. 



because the rescattering plays a role that it decreases i>2 
and tends to make the reaction system isotropic from 
above, which also is consistent with the isotropic fea- 
ture of rescattering angle in LUCIAE rescattering model. 
So the anisotropy of momentum space, i.e. v^, may be 
produced in partonic state before hadronization if the 
hadronic rescattering has no direction-sense. In other 
words, flow must be formed before the hadronization. 

(Ill) We only remain initial anisotropy of momentum 
space but lack initial anisotropy of coordinate space when 
hadronization takes place as shown with triangles in fig- 
ure 5, the rescattering still decrease hadrons' V2 and creat 
a big anisotropy of coordinate space €2- 

Not only in above three conditions but also in our de- 
fault condition, we all find that hadrons' rescattering pro- 
duces many secondary hadrons and suppresses the ellip- 
tic flow V2 of hadrons to a certain extent. When the two 
anisotropy of coordinate and momentum space tend to 
be homologic with time, the rescattering ends and the 



reaction system freezes out. 

However, we note that the effect of rescattering on V2 
in high px domain is not perfect from Figure 3, i.e., 
V2 does not overlap with experimental data which sat- 
urates above Pt ~2 GeV/c. We think that one cause 
may be from our anisotropic factor in the momentum 
space. Even though the experimental data can be fitted 
well in lower px, it is a bit artificial since the nucleon 
distribution of transverse plane (x,y) and anisotropic hy- 
drodynamical evolution are not taken into account in this 
hydrodynamics The other cause could be the par- 
tonic effect on flow, such as partons cascade [2f| and jet 
quenching |27l |28| , is also very important for the charac- 
ter of particles with high px, which we also do not take 
into account. Of course, the main goal of this work is to 
investigate the hadronic rescattering effect on the final 
state hadronic elliptic flow. 



CONCLUSIONS 



In conclusion, we apply hydrodynamics model followed 
by the rescattering model to investigate the effect of 
hadronic rescattering on elliptic flow v 2 for 200 GeV/c 
Au + Au at 20-40% centrality. Hadronization is treated 
by Cooper-Frye method and an initial anisotropic phase 
space of hadrons is assumed. We find that the yield of 
hadrons is increased after the rescattering. While, the 
rescattering among the hadrons plays a suppression role 
for t>2 , which makes an asymmetric system in momentum 
space less anisotropic. The suppression becomes weaker 
with the increasing of transverse momentum. We miss at 
least 20 ~ 40% V2 after hadronic rescattering in our hy- 
drodynamics + rescattering model. This may, however, 
change if hadronization occurs when partons make their 
last scattering as treated in Texas multi-phase tr ansp ort 
model when the spatial anisotropic is small [2j| l30j in- 
stead of when the energy reaches certain critical energy 
density presented here. On the other hand, the hadronic 
rescattering makes the coordinate space of hadrons tend 
to the similar directions to the momentum space. In our 
model, we also find that the rescattering does not induce 
the elliptic flow if only initial anisotropy of coordinate 
space exists in the hadronic phase. 

Considering that the hadronic system created in RHIC 
collision is very dense, the following hadronic rescatter- 
ing may happen potentially and play an important role. 
Of course, only information of the momentum space af- 
ter the freeze-out can be measured in experiment. But 
anyway, the effect of elliptic flow V2 from final hadronic 
rescattering deserve attention if the hadronization takes 
place when the energy reaches certain critical energy den- 
sity presented here in hydrodynamics. 
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